Inferring effective interactions from the local density of states: application to STM 

data from Bi2Sr2CaCu208+5 
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While the influence of impurities on the local density of states (LDOS) in a metal is notoriously 
non-local due to interference effects, low order moments of the LDOS in general can be shown to 
depend only on the local structure of the Hamiltonian. Speciflcally, we show that an analysis of 
the spatial variations of these moments permits one to "work backwards" from scanning tunneling 
microscopy (STM) data to infer the local structure of the underlying effective Hamiltonian. Applying 
this analysis to STM data from the high temperature superconductor, Bi2Sr2CaCu20g+5, we find 
that the variations of the electro-chemical potential are remarkably small {i.e., the disorder is, in a 
sense, weak) but that there are large variations in the local magnitude of the d-wave gap parameter. 
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I. INTRODUCTION 

Scanning tunneling microscope (STM) measure- 
ments on highly correlated materials, such as the 
high temperature superconductors, are increasingly be- 
ing gnal'-zc d^-^-^-^-^-^-^-^-^-^"--^^-^^--^^-^'^--^^-^^--^^-^^'-^^-^"-^^-^^ 
in the hopes of inferring information concerning the char- 
acter of the dominant interactions affecting the electron 
dynamics. However, it has long been known, as exem- 
plified by a series of famous experimentsSSiSi on simple 
metal surfaces, that the effects on the local density of 
states (LDOS) of a change in the Hamiltonian are highly 
non-local. The interference patterns produced by quasi- 
particles scattering from an impurity cause ripples in the 
LDOS to spread far and near. The possibility of such 
non-local effects make it a non-trivial issue how to best 
solve the inverse problem: given a set of high quality 
STM data, how does one work backwards to infer the 
character of the effective Hamiltonian that produced the 
measured spectra? In the present paper, we show that 
the low order moments of the LDOS are determined en- 
tirely by the local form of the Hamiltonian, and that they 
can therefore be simply analyzed to obtain an approxi- 
mate solution to this inverse problem. 

Of particular interest to us are STM studies of 
the high temperature superconductor, Bi2Sr2CaCu208+5 
(BSCCO), which have revealed^i^^i^i^i^dOi^d^d^ intrigu- 
ing spatial inhomogeneities in the LDOS, Pr{(^)- For en- 
ergies to of order the superconducting gap, Aq, there are 
order one variations in the LDOS, although for \uj\ ^ Aq, 
and again for |aj| ^ Aq there are only relatively small 
fractional variations in p as a function of position, r. The 
energy of the low energy peak in Pr(w) (usually called 
the superconducting coherence peak) varies substantially 
from one location to another, with a "patch" diameter of 
around 30A. There have been at least two different inter- 
pretations of these results proposed: The first proposed 
interpretation^iH associates the local value of the peak 
in Prii-o) with the local magnitude of the gap parameter, 
Ar.r', in a mean field Hamiltonian. However, Fang et al.^ 
showed that this interpretation can be quite misleading. 
They calculated the LDOS of a metallic region (with a 



local value of A = 0) embedded in the background of 
a d-wave superconductor, with a bulk A = Aq. Unsur- 
prisingly, they found that for a patch size comparable to 
the superconducting coherence length, the proximity ef- 
fect causes the LDOS in the metallic region to look like 
that of a uniform d-wave superconductor, with peaks at 
nonzero energy. (See also Ref. 5.) 



II. LOW ORDER MOMENTS OF THE LDOS 

Moments of the one particle spectral function can be 
expressed in terms of matrix elements of the Hamiltonian 
and certain thermodynamic correlation functions, as has 
been recognized at least since the work of Nolting et alm^. 
Nolting et al. and later authors, including Oganesyan^— , 
Randeria et al^'^ and Norman et al."^^ proposed using a 
moment analysis to extract information from angle re- 
solved photoemission (ARPES) data. While the moment 
analysis can be carried through at a fundamental level, 
i.e. for Schrodinger's equation, we are actually more in- 
terested in interpreting the results in terms of a low en- 
ergy effective theory. We therefore begin by deriving ex- 
pressions for the first and second moments of the LDOS 
in the presence of disorder (which was not considered by 
Nolting et at), and for different assumptions concerning 
the structure of the effective Hamiltonian. 34 

We define the moments of the LDOS Pr{^) as 



M„(r) = 



(2.1) 



(Note that an STM experiment does not actually mea- 
sure the LDOS due to the existence of unknown matrix 
elements. However, since we have normalized the mo- 
ments as indicated, the effect of any energy independent 
matrix element cancels in M„.) By definition, Afo(r) = 1. 
Expressions for the first and second moments in the two 
cases discussed here are derived in Appendix A. 

Case I: A quadratic (mean-field) effective 
Hamiltonian: In systems in which the relevant physics 
can be described in mean-field approximation, the low 
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energy physics can be derived from an effective Hamilto- 
nian which is quadratic in the fermionic field operators; 
this is the case in which the moment analysis is most pow- 
erful, as the moments depend only on the Hamiltonian, 
and are independent of the associated thermodynamic 
correlations. 

Needless to say, this approach neglects much of the 
strong correlation physics. We partially address this issue 
in Case II below. For now let us consider the most general 
mean-field Hamiltonian, 

H = -\Y. [Uy^Hr' + ^r,M4' + H-c.} , (2.2) 



for which the moments are straightforwardly seen to be 

Ml(r) = -<r,r = -Mr (2.3) 



M2(r) = ^Wir',r + Xl 



(2.4) 



Here the index r labels the position and band index of 
a Wannier function, as well as the spin-polarization, /ir 
is the local value of the electro-chemical potential, t^y 
is the hopping parameter from r to r', and Ar.r' is the 
superconducting gap parameter. 

Case II: Coupling to fluctuating fields: In Case I, 
we assumed that A, /i and t were simply numbers. How- 
ever, these parameters could have dramatically fluctuat- 
ing pieces, which are either present due to phonon medi- 
ated deformations of the effective electronic Hamiltonian 
or are generated (via a Hubbard-Stratanovich transfor- 
mation) from an underlying interacting Hamiltonian. As 
discussed at the end of Appendix A, all the expressions in 
Eq. (|2.3|l and (|2.4|) are generalized in this case by replac- 
ing the left-hand side of each expression by its thermal 
expectation value (signified by (••■)), so in Eq. 1)2. 3|l . 
/ir — > (/ir) and in Eq. (|2.4|l . ir,r'ir',r — > {tr,r'tr',r) and 
Ar,r'P (Ar,r'A]!, J.). For example, a fluctuating su- 
perconducting gap parameter with rms value |A| (above 
Tc) makes the same contribution to the second moment 
as a mean-field gap of magnitude |A| (in the supercon- 
ducting state). 

Henceforth, we will focus on Case I, with the under- 
standing that the effective parameters that enter the 
Hamiltonian may inherit temperature and magnetic field 
dependences from the underlying thermodynamic state 
of the system. 



III. PARTIAL MOMENTS 

The integrals which define the moments, M„(r) extend 
to infinite energies. Not only does this make them un- 
measurable quantities, even if we could measure them, 
we would not expect them to yield useful information. 
An effective Hamiltonian provides a reasonable descrip- 
tion of the complex microscopic problem, at best, over a 



finite range of "low" energies. Given this, we need to con- 
sider the more complex issue of what information can be 
gleaned from the partial moments, in which the integrals 
are cutoff at a finite energy f2: 



M„(r;r!) 



Ul"pr{uj)dcj 



(3.1) 
(3.2) 



Note that M„(r) = M„(r;oo). 

There is no simple general expression for the partial 
moments. However, there are circumstances in which im- 
portant information can be obtained from them, anyway. 
In this sense, the application of these ideas is analogous 
to the application of the /-sum rule to the analysis of 
optical conductivity data. For instance, were we inter- 
ested in the temperature (T) or the magnetic field (B) 
dependence of the effective Hamiltonian, we might look 
at the difference between the LDOS at T = and B = 
and its value at non-zero T or B; if, as is often the case, 
Pr{ui) is essentially T and B independent above a charac- 
teristic energy, 51, then by analyzing the partial moments 
M„(r, f2), we can obtain information about the induced 
changes in the effective Hamiltonian. 

For the present purposes, we will focus on a different 
example, relevant to the analysis of the above mentioned 
STM data in BSCCO. Suppose we are interested in ob- 
taining information about the spatial variations of the 
effective Hamiltonian. Again, assuming that for a range 
of energies uj > Q, Pr{i^) is approximately r independent, 
we can analyze the spatially varying part of the moments, 
5Mn = Mn{v,VL) - Mn{^) whcrc M„(0) is the average 
of M„(r, rj) for all r, in terms of the spatially varying 
pieces of the effective Hamiltonian. An added advantage 
of this approach, as we shall shortly see, is that differ- 
ent terms in the effective Hamiltonian produce changes 
in the moments on different energy scales. 

An explicit model: To illustrate the various aspects 
of this approach, we will consider a specific model prob- 
lem, which we design with the STM studies of BSCCO 
in mind. We suppose that the electronic structure is gov- 
erned by a 2D effective Hamiltonian, H'^-^-f on a square 
lattice of the form given in Eq. H2.2|l . with parameters 
taken to be representativeiSiSS of the band structure and 
pairing symmetry of BSCCO. We consider the case in 
which there is an inhomogeneous "patch" of size 6x6 
embedded in an otherwise uniform infinite background. 
(While this choice of patch size is arbitrary, we were mo- 
tivated by the experimental observation in BSCCOi^ifiiii 
of apparent patches of size around 3nm.) For the spa- 
tially uniform piece of H'^^^ we define ir,r' = t for r 
and r' nearest-neighbor sites. (The band structure re- 
tains up to fifth nearest neighbor terms, with second 
nearest neighbor ir.r' = — 0.55i, and other hopping ma- 
trix elements smaller in magnitude by a factor of five 
or more.) The chemical potential is tr,r = /-* = 0.44i, 
and Ar.r' = iAg for r and r' nearest-neighbor sites with 
the signs chosen consistent with the d-wave symmetry of 
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the order parameter. We further take a representative 
vahie of Aq = 0.085 i. (In the case of BSCCO, t « 0.15 
eV and Aq « 0.04 eV.) The band, which is not particle 
hole-symmetric, runs from —2.7t to 6.7t. 

Inside the patch we will consider three different ex- 
treme cases: (1) Inhomogeneous electro-chemical poten- 
tial: /ir = ^ + Sn for r inside the patch, (2) Inho- 
mogeneous gap parameter with the d— wave symmetry: 
Ar,r' = Aq ± SA ioi nearest neighbor r, r' inside the 
patch. 3) Inhomogeneous band structure: ir,r' = t + St 
for nearest neighbor r, r' inside the patch. 



-0.4 

to 

-o.s 



^ -1.2 

to 






6.7 ( 




3.4 ( 




1.4 ( 




0.34 r 



45M,(i- =00) 



region of suppressed A 

^^^^^^^^^^^^^ 




48M ir = 0) 



position r 

FIG. 1: (Color online.) The dimensionless response 
SM2{r; n)/5A^ due to the suppression of A in a 6 x 6 patch 
(displayed region is 14 sites wide) for several choices of fl. 
(The band runs from —2.7t to &.7t.) As we expect, when the 
cutoff is large, the second moment only knows about the local 
environment, and its value changes abruptly as r crosses into 
the region of suppressed gap. Surprisingly, the response of 
M2(r;f2) to spatial variations of A remains quite local for Q, 
well above Ao = 0.085t. 

An exact formal expression for the single-particle 
Green function of this problem can be obtained from 
multiple scattering theory in terms of the Green func- 
tion of the uniform system. For each of the three cases, 
we have evaluated these expressions numerically to com- 
pute the LDOS. These are in turn integrated numerically 
to compute the partial moments. 

Locality of the partial moments: Fig. ^ shows 
the response of M2{y]VI) to the inhomogeneity in the 
gap parameter for several choices of il. In the present 
case, we have taken 5/S. = — Aq, that is we have con- 
sidered the case in which the pairing potential is absent 
inside the patch. For presentational clarity, we have plot- 
ted 4(5M2(r; ri)/5A^, the renormalized difference between 
M2(r; ri) and the average value M2(ri), which in this case 
is identical to the asymptotic value lim|r|_»oo -^2(r; ri)). 
Hence 5M2{v; Vl) approaches when the position r is far 
from the patch. Both the LDOS at fixed lo (not shown) 
and M2(ri) for O < Aq exhibit complicated r depen- 
dences, with induced spatial variability that extends far 
outside the patch. However, once is greater than a few 
times Aq, even though Q, is still substantially less than 
the bandwidth, M2(r; fi) approaches its asymptotic value 



it follows that, in the limit — > cxd, 
^ -0.18^2. For non-infinite VL, 



rapidly. 

From Eq. 

<5M2^-^- — ^ 
we can see that the partial moment is already approxi- 
mately equal to the 17 — > cxd moment when w 3i (i.e., 
half the bandwidth) and it overestimates the response of 
the moment by w 40% for Q. = lAt {i.e., 1/5 of the band- 
width.) Thus, qualitatively accurate information about 
the structure of the effective Hamiltonian can be obtained 
even if the range of integration is restricted to as little as 
1/5 of the bandwidth. Specifically, if one were to simply 
pretend that 5M2{v\ Q. = 1.4i) is approximately equal to 
(5M2(r; = cxd), one would conclude that the gap param- 
eter changes from its value inside the patch to its value 
outside the patch over an interface region of width equal 
to a couple of lattice constants (as opposed to the actual 
one lattice constant) and one would underestimate the 
magnitude of the variation of the gap parameter by 40%. 



Disti ngui shi ng in homogeneities in A, /i and t: 

and 



In Figs. 2(a) and 2(b) we have plotted the spatially vary- 
ing part of the first and second moments, 6Mj{Y,Vl) = 
Mj(r, il) — Mj{Q) for each of the three different types 
of inhomogeneity as a function of the cutoff, Vl. Here, r 
is the point at the center of the patch. Specifically, the 
three curves are for the case of inhomogeneous electro- 
chemical potential (dotted lines) with 5^ = 0.1 /x, inho- 
mogeneous gap parameter (solid lines) with (5 A = — Aq 
(i.e. vanishing gap parameter inside the patch) and in- 
homogeneous band structure (dashed lines) with 5t — 
0.1 1. The curves have all been normalized, as indicated, 
to be appropriate dimensionless measures of the varia- 
tion. Since the spectrum is bounded above, the moments 
achieve their asymptotic values once $7 is in excess of the 
band edge, 6.7t. The various kinks in the curves reflect 
the lower band edge. 



In Fig. 2(a) we see that in the case of an inhomoge- 
neous /i(r), the first moment has already reached 40% of 
its asymptotic value by the time is 1/20 of the band- 
width. The effect of the inhomogeneities in A and t goes 
to zero for O in excess of the bandwidth as we expect 
from Ea. H2.3|l . This happens very rapidly for the case 
of gap inhomogeneities, so even when f2 is moderately 
small compared to the bandwidth, the first moment gives 
a reasonable measure of the inhomogeneity of the band 
structure parameters, and t. 



The most striking feature of Fig. 2(b) is that for small 
ri, 5M2{v\ il), becomes rapidly large (even larger than its 
asymptotic value) in the case of an inhomogeneous A, but 
remains relatively small in the other two cases until is a 
substantial fraction of the bandwidth. (This is especially 
true for inhomogeneous fj,.) We therefore conclude that 
(^M2(r, ri) for Vl a few times Ag but still much less than 
the bandwidth gives a reasonable measure of the spatial 
variations of A. 
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1 Sh- = 6(f) 
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FIG. 2: (Color online.) The diflterence between partial mo- 
ments inside and outside the patch in response to the inho- 
mogeneity in the parameter of the Hamiltonian h — A, as 
a function of the cutoff il. Sh = 10 meV in all figures, (a) The 
first moment difference (Eq. 12.31 '). (b) The second moment 
difference (Eq. lEU). 



A. LDOS near a unitary scatterer 

We have also analyzed an example of a single strong 
(unitary) scatterer imbedded within an otherwise homo- 
geneous superconductor. The existence of a resonance 
near zero energy in the presence of a unitary scatterer 
has been predicted by Balatsky et al.^ and experimen- 
tally confirmedi^ii^*^. Also finite energy resonances in 
the presence of a magnetic (non-unitary) impurity have 
been observed in Refs.^"*^^ 

Fig.lHJa) displays the LDOS p{uj) at different positions 
away from a unitary scatterer along the lattice direction 
(gap antinodal direction). There is subtle fine structure 
in the LDOS which we have suppressed by introducing an 
energy broadening 6u! =2 meV. As is clear in the Fig.|2Ia- 
b), there is a w « resonance that peaks one site away 



FIG. 3: (Color online.) (a) LDOS p{ijj) near a unitary scatter- 
ing impurity for the model introduced in Sec. Illll at different 
positions, (b) Log of the peak height at « as a function 
of distance from the impurity, (c) SMi (r; f2) along the bond 
direction for several choices of fl. 



from the impurity (pr=o(^^ = 0) = 0); ^-^^d decays over the 
span of several sites. In the direction 45° to the nearest- 
neighbor bonds (nodal direction), the peak intensity is 
substantially lower at short distances (see Fig. 2)). As is 
seen in the Fig. ISJc) , the first (partial) moment is much 
more highly localized, even at relatively small Q, despite 
the fact that the impurity site LDOS weight has been 
pushed almost entirely beyond the integration range. 

As an aside, we note that it may be surprising that 
the disturbance caused by the impurity appears to ex- 
tend more strongly in the antinodal direction than in 
the nodal direction, whereas it is theoretically clear that 
the asymptotic effect of the impurity at long distances is 
stronger in the nodal than in the bond direction. How- 
ever, explicit calculations reveal that, to a remarkable 
extent, the structure at short and intermediate distances 
depends strongly on the band structure. This sensitivity 
to band structure is illustrated in Fig.^ We have plotted 
the LDOS around a unitary scatterer, for a band struc- 
ture e(k) = erm(k) -I- Q!(eo(k) — e„„(k)), which interpo- 
lates between a simple nearest-neighbor band structure, 
e„„(k) (at a = 0), and the full band structure, eo(k) (at 
a = 1), used throughout this article.^^ As a increases 
from to 1, the short range character changes from pre- 
dominantly nodal to antinodal falloff at short distances. 
Of course, the falloff is always slower in the nodal direc- 
tion at long distances for any a £ [0, 1]. 

We have also looked at the case in which a resonant 
scatterer is placed at the center of a patch with either 
reduced or enhanced gap. The behavior of the impurity 
induced disturbance of the LDOS is not qualitatively dif- 
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FIG. 4: LDOS at = 0, near a unitary scatterer for several 
band structures, e(k) — e„„(k)+Q(eo(k) — enn(k)) on a region 
of size 30 X 30 sites. [Lighter shades denote larger p{r, ui — 0).] 
Here, a = corresponds to nearest- neighbor tight-binding, 
whereas a = 1 corresponds to the full band structure from 
Refs. [1^12911 used throughout this article. The short-range 
character of p(r, = 0) is sensitive to band structure, whereas 
the long-distance disturbance is strongest along, or near, the 
nodal direction. 



as indicative of the existence of regions with different 
superconducting character, is largely reproduced in the 
second moment map, but seems largely uncorrelated with 
features in the topograph or first moment map. 




(a) (b) 

FIG. 5; (Color online.) (a) The topograph of BSCCO sam- 
ple, (b) The first moment map A/i (r; 200 meV) of the same 
sample. (Units in meV.) The lattice superstructure is clearly 
visible in both. 



ferent in these cases, either - there is still a resonant peak 
near uj = and its spatial decay occurs with a similar 
pattern. Of course, other features of the LDOS are even 
more complicated in these cases. Nonetheless, even not- 
ing the strongly non-universal character of the induced 
features in the LDOS at short distances, the fact that the 
low energy peak observed in STM studies of Zn doped 
BSCCO is sometimes seen to fall by as much as a factor 
of 100 in intensity within two or three lattice constants of 
a Zn impurity-'^'*'^^ is not easily reconciled with a simple 
quasiparticle theory. 



IV. BSCCO ANALYSIS 

As a test of the practical usefulness of the local moment 
analysis, we have applied it to some STM data of Fang 
and Kapitulniki on a 220 x 280A section of surface of near 
optimally doped Bi2Sr2CaCu208+5. Fig.Elshows a plot 
of (a) the topograph and (b) the first moment Mi(r; fi) 
with a cutoff fl =200 meV. (We chose this cutoff because 
it is several times the average gap maximum but still 
small compared to the bandwidth, and because at this 
energy the LDOS is more or less spatially uniform as 
Fig. [3 shows). Fig.Elshows a plot of (a) the gap map and 
(b) the second moment also with — 200 nieV. (The gap 
map displays the energy difference between the positive 
and the negative energy maxima in the LDOS at each 
point r.) The superstructure due to the buckling of the 
Bi-0 planes is clearly evident in the topograph and the 
plot of Ml (Fig. [SJ but not in the gap map or the plot 
of the second moment (Fig. (HI. Conversely, the splotchy 
character of the gap map, which has been interpreted 




FIG. 6: (Color online.) (a) The gap map of BSSCO 
near optimal doping, (b) The map of the spatially vary- 
in g piece of the gap mag nitude from the second moment, 
y'Ahir; 200 meV) - A'/j"'" from the same sample. (In deter- 
mining M™*", the spatial minimum value of A/2(r), we have 
ignored values of Af2(r) that are smaller than the mean by 
more than three standard deviation.) 

On the basis of the analysis of the previous section, 
we expect that Mi(r, O) is approximately proportional 
to the electro-chemical potential, even though fl is small 
compared to the band width. One of the most notable 
conclusions is that Mi varies only by about 20 nieV from 
position to position. Moreover, most of that variation is 
due to the superstructure variation. Even if we take into 
account the fact that having a cutoff of 200 meV may 
mean that the spatial variation in the electro-chemical 
potential, Sfj,, may be as much as twice the variation in 
Ml, we still conclude that, in units of the band width 
(1.4 eV), the variations in the potential are small. 

The second moment varies as a function of r by about 
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7: The spatial variance of the normahzed LDOS, 
is calculated and plotted as a measure 



p(") 



of degree of inhomogeneity of LDOS energy range between 
— n = —200 meV and il = 200 meV. As you can see the normal- 
ized LDOS is very homogeneous outside of the energy range 
between — lOOmeV and 100 meV. Also, as has been noted be- 
fore, the LDOS is homogeneous at low energies as you can see 
in the flat region between — 25meV and 25meV and the sur- 
prisingly sharp onset of the inhomogeneity at about ±25 meV. 



3000 (meV) 2. On the basis of the results in Fig. |2(a) 
we feel it is safe to associate this with variations of |A|^ 
However, because of the relatively small value of fl, direct 
application of Ea. l2.4l mav overestimate the variations of 
Ap by as much as a factor of 2. 

Assuming this factor of 2 overestimation i.e., 
(5M2(r;200meT/) « 2SM2{r) w 2^ and that the re- 
gions with the smallest values of M2(r, $7) correspond to 
places where A(r) = 0, this analysis leads to the con- 
clusion that A varies as a function of position from 
to about 80meV; (Fang et ali- have shown the experi- 
mental data is compatible with minimum A = 0.) Were 
we to assume a non-zero value for the minimum value of 
A, we would obtain a correspondingly larger estimate of 
the maximum gap. In comparison, were we to interpret 
the gap map as representative of A(r), the range of val- 
ues inferred would be from approximately 25 to 70 meV. 
While both methods of analysis suggest large amplitude 
variations in A, we believe that the present estimate, 
which implies even larger variations than does the gap 
map, is the more reliable. This exaggerated response of 
the electronic correlations to small variations in the local 
potential is suggestive of some form of frustrated phase 
separations^. 
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V. APPENDIX A 

Here we outline the proof for the sum rules, first for the 
mean field case. The mean field Hamiltonian of Eq. 12.21 
can be written as 



H — £(73 -I- AtTl, 



(5.1) 



in which e and A are in general nontranslationally in- 
variant hopping and gap matrices, and aj are the Pauli 
matrices with Nambu indices. The Green function for 
this Hamiltonian can be defined as 



G = 



1 



(5.2) 



{uj + i6)I - H' 
In terms of G, the local density of states can be found as 



Pr{uj) ^ --Im [Gi,i(r,r)] 



(5.3) 



where the subscripts are the Nambu indices. We can 
decompose G into the sum of a Hermitian and an anti- 
Hermitian operator. The Hermitian part obviously does 
not contribute to the LDOS. The anti-Hermitian part is 
equal to 



1 



(G-G^) = -i7rS{ujI - H). 



(5.4) 



Of course S{ujI — H) is defined as a matrix with the same 
eigenvectors |A) as H, but with eigenvalues 5{ljj — A). 
Eq.O and Eil tell us 



From this, we get the following relations: 



However, 



(5.5) 



TTtl 

— '"rl;rl- 

^Tl;rl = lrl;rl = 1 
^rl;rl = erl;rl = ^ M 



2 

rr' 5 



(5.6) 

(5.7) 
(5.8) 
(5.9) 

(5.10) 
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and this gives us the desired sum rules. 

For the general interacting Hamiltonian, we can get the 
density of states from the imaginary part of the retarded 
Green's function: 







Im 




TT 




1 






e 


2n 7_oo 



^^'{{i,l{0),Mt)})dt, (5.11) 



and therefore 

oo 

aj"^ Pr{Lu)du; 

j_ r°° \f_.£ 
1 

2^ 



6{t) 



{{Mt),ijl{0)})dtdu; 



from which the most general counterparts of Eq. (|2.3() and 
(|2.4|) follow. As an application of this formula we can 
see that if in the Hamiltonian 1)2. 2|l . t^.r' and Ar.r' are 
turned into dynamical fields by adding a term to the 
Hamiltonian that is a function of these fields and their 
momenta then Eq. (|2.3(l should be replaced by 



Mi(r) = -(A^r) 



(5.13) 



(5.12) 



in which {fi^) is the thermal average of the dynamical 
field fir- For the second moment, Ea. (|2.4fl turns into 



M2(r) = ^(tr,r'ir',r)+^(Ar,r'A]:,_^) 

J./ J,/ 

+ {[tr,r,H]). (5.14) 



The expectation value of the commutator of the Hamilto- 
nian with any operator in equilibrium is, of course, zero, 
so the last term drops out of Eg. H5. 14(1 . 
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aZi— , based on a "Mottness" map which, within a single- 
band Hubbard model, relates certain integrals of the LDOS 
to the expectation value of the local electronic density, 
n(r). 



The nearest-neighbor hopping in e„„ (k) has been set larger 
than in eo(k) to insure the bandwidths are the same for 
both a = and a = 1. 



